Transient effects in the backscattered current of a Luttinger liquid 
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We study the backscattered current in a Luttinger liquid in the presence of a point like weak 
impurity switched on at finite time, taking into account finite-temperature effects. We show how 
the well-known results for a static impurity are distorted. We derive a dimensionless parameter tr 
as function of the electron-electron interaction and the temperature, such that for tr < 1 (> 1) 
the switching process is relevant (irrelevant). Our results suggest the possibility of determining the 
value of the Luttinger parameter K through time measurements in transport experiments at fixed 
voltage. 
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Ph ; I. INTRODUCTION 

\o 

Quantum transport in novel one-dimensional (ID) materials, such as quantum wires and carbon nanotubes^, is one 
of the most active areas of present research in condensed matter physics. Recent experiments have confirmed some 
of the striking effects that characterize the Luttinger liquid (LL) picture of ID nanostructures 2 -, such as spin-charge 
^ . separation 3 - and charge fractionalizatiorA One central issue in all experiments is the determination of Luttinger 
exponents that depend on the rigidity constant K, which parametrizes electron-electron (e-e) interactions. These 
exponents have been measured in pioneering tunneling experiments 5 . However, some doubts still remain in the 
> ' determination of K due to the fact that dynamical Coulomb blockade also leads to power laws in systems surrounded 
by ohmic devices^. Moreover, it has been shown that environmental resistance can also contribute to the measured 
value of K-. It is then very important to conceive alternative ways to measure K . To a large extent, quantum transport 
| is highly non trivial due to the influence of a variety of combined effects produced by impurities, junctions, contacts, 
etc. In particular, in the study of impurities in Luttinger liquids, an observable of special interest is the backscattered 
current lbs- For a point like static impurity, lbs opposes the background current and goes as V 2K ~ 1 , where V is 
the bias voltage^. The effects of finite temperature and finite length of the quantum wire lead to characteristic non- 
monotonic parameter dependencies of the current and yield a rich structure in the noise spectrum which depends on 
I \ . For a time-dependent oscillatory impurity, the current grows and the conductance of a one-channel quantum wire 

■ is greater than its background value e 2 /h for strong repulsive interaction (K < 1/2) 10 11 . More recently, the effect of 
t-H | several time-dependent impurities was considered at zero temperature and infinite lengt h 12 13 . When temperature is 
taken into account, for two oscillating barriers, the previous zero temperature results are distorted; the behavior of 
the backscattered current and the change in the differential conductance depend on different regimes which can be 
established as function of the frequency, the temperature and the separation between the impurity potentials^. 
' One aspect that has been seldom explored in the context of the Luttinger liquid is its response to a sudden switch of 
an interaction of the system with an external field. As an example, the total energy density of a Tomonaga-Luttinger 
liquid in the presence of a forward-scattering time-dependent impurity switched on at finite time has been computed 
exactly, distinguishing two well defined regimes in terms of the relationship between the frequency of the perturbation 
and the electron energy. This study allowed to determine a time interval in which the finite-time switching process is 
k^j ' relevant. 

In this paper we examine the effect of a transient process on transport properties of a Luttinger liquid. When a 
local barrier, that can be considered as a backscattering impurity, is turned on at a finite time to, a backscattered 
current is produced. Let us stress that such a barrier can be created, for instance, in a single- walled carbon nanotube 
by applying a voltage to a narrow metal gate electrode^ (See Fig. I). We obtain an analytical expression for Ib s as 
function of time. The envelope of this function decays in a way that crucially depends on e-e interactions, i. e. on 
the K parameter. This allows us to find a simple relation between the decay process and the value of K. Thus, we 
show that K could be determined by measuring time intervals within the reach of recently developed pump-probe 
techniques with femtosecond-attosecond time resolutions^. Then, in contrast to the conventional techniques used 
up to now to measure K , based on the study of stationary transport properties for varying voltages, we indicate 
an alternative way based on the analysis of the transient current for fixed voltage. Apart from its interest in the 
context of strongly correlated electronic systems, our studies could be relevant in the area of cold atomic gases, where 
quantum quenches are being intensively investigated^. 

The work is organized as follows. In Section II we present the model and review the calculation of the backscattered 
current in terms of a vacuum expectation value of exponentials of bosonic fields. In Sections III and IV we present the 
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FIG. 1: The figure shows a quantum wire coupled adiabatically to two reservoirs with different chemical potentials, with a 
backscattering impurity switched on at time t — to. The current I(t) is measured as a function of time. 



results obtained at zero and finite temperature, respectively. We define a dimensionless relaxation parameter tr as 
function of the temperature and of the electron-electron interaction. Finally, we analyzed and summarize our results 
in Sec V. 



II. THE MODEL 



We consider a long and clean LL adiabatically coupled to two electrodes with different chemical potentials, such 
that fJ-L — fJ-R = eV. We will restrict our study to the case in which the electrodes are held at the same temperature. 
This condition is very important in order to apply standard bosonization techniques^. Indeed, as was very recently 
explained 20 , standard bosonization is expected to work well only in special situations corresponding to small deviations 
from equilibrium. The non trivial ingredient of our model is a barrier of constant height (an externally controlled 
impurity) that is turned on at a finite time to, producing backscattering of incident waves. Under these conditions 
we model the LL in terms of the following Lagrangian densitjJ^: 



L — Lq + Li mp ) (1) 

where 



1 / d 2 d 2 \ 

L Q = -^(x,t)[v 2 --^J^x,t) (2) 

describes a spinless Tomonaga-Luttinger liquid with renormalized velocity v and 

Limp = — ^tHx - x )Q{t - t Q ) cos \2k F xlh + 2VttKv^(x, t) + eVt/h] (3) 
7mA L 

represents the scattering of spinless electrons with the external barrier at the point xq and switched on at the time 
to, with a coupling amplitude gi>. V is the external voltage applied to the quantum wire and K measures the strength 
of the electron-electron interactions. For repulsive interactions K < 1, and for noninteracting electrons K = 1. A 
is a short-distance cutoff. In the above expression we only take into account backscattering between electrons and 
impurities, because the forward scattering does not change the transport properties studied here, at least, in the 
lowest-order of the perturbative expansion in the couplings. 

In the absence of the impurity, the background current is Iq = e 2 V/h . When the impurity is acting the total 
current is I = Iq — lbs ■ The operator associated to the backscattered current is defined as^ 



hs(t) = ^rr®(t - t ) sin[2k F x a /h + 2VirKv$(x , t) + eVt/h] . (4) 
7rnA 



The backscattered current at any time t is given by 



I bs (t) = (0\S(-oo;t)T ba {t)S{t;-oo)\0), 



(5) 
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where (0| denotes the initial state and S is the scattering matrix, which to the lowest order in the coupling gs is given 
by 



S(t; — oo) = l — i 



dx 



(6) 



When one inserts ((6]) in (|5]) one finds several terms of the form 



A a ^ = (0| exp[2iaVTrKv<f>(x, t')] exp[-2if3VirKv<£>(x, t)]\0) , 



(7) 



with a, (3 = ±1. This kind of vacuum expectation values (v.e.v.) of vertex operators have been computed many times 
in the literature. It is well-known that A a ^ a = 0. Thus, the building block of our computation is A a ^ a . Let us sketch 
the calculational procedure for a = 1. Using Baker-Haussdorff formula and taking into account that the commutator 
of the fields is a c-number, we can write: 



A M = (0| exp [2WitKv{<S>(x, t') - $(a;, t))] |0) exp [2ttKv[^(x, t'),$(x, t)]] . 



(8) 



At this point we observe that the v.e.v. of the exponential in the first factor above, can be written as the exponential 
of a v.e.v.: 



(0| exp [2W-kKv($(x, t') - 3>(x, t))] |0) = exp [ - 2irKv (0|($(x, t') - t)f |0)] . (9) 
Now, in order to explicitly evaluate the previous expressions we need Keldysh-l lesser function G < given by 



(0|$(M')$OM)))|0) =iG<(x,t;x,t') = A- J dpdtoer^- 1 '^ (8(-uj) + n B (\oj\)) 5(ui 2 - v 2 p 2 ), (10) 
where nB(|w|) = (e k B T — 1)~ • Putting all this together we finally obtain 



(0| exp^vVftTu^x, t')] e^p[-i2VTrKv^(x 7 1)]0) 
where sgn is the Sign function. 



(Kirk b T/hv) 2K exp[iTrKsgn[t - t'}} 



\smh[Trk b T(t-t')/h]\ 



2K 



(11) 



III. RESULTS AT ZERO TEMPERATURE 



First of all, we compute (|5|) at zero temperature and with the impurity switched on at time — oo 



I bs (t, -oo) 



<?|eA 



2K-2 



2irh 2 v 2K T[2K} 



eV 

h 



2K-1 



sgn[V]. 



(12) 



This corresponds to the well-known case of a static impurity^, where the backscattered current goes as V 2K ~ 1 . We 
note that for K < 1/2, the backscattered current becomes large when V decreases. Hence, the perturbative expansion 
in powers of gs breaks down when V — > 0. Using a scaling analysis we can estimate that this expansion is valid when 
j^;{^^-) K ^ 1 *C 1. We emphasize that expression (TT21 does not include the case V = 0, where the current is zero 
too. All these statements imply that the current must be a nonmonotonic function of V. In order to determine this 
function one has to go beyond the lowest-order perturbative results of this work . 

Now, we compute the backscattered current when the impurity is switched on at a finite time to: 



hs{t,to) — 



G(t)T[2K\ 



-2-2K 



T[K]T[2 - K] 



iF 2 [l - K; 3/2, 2 - K; -(r/2) 2 ] I bs (t, -oo) 



(13) 
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FIG. 2: Example of an experimental plot of / as function t, leading to the determination of K. We have set = 10 12 (l/seg) 

(corresponding to V ~ mV), V- = (O.l) 1 ^ 2 , = 1. I is measured in units of ^-rf-- For this particular case, using the envelope 
(dashed line) one obtains K = 0.4. 

where iF% is the generalized hypergeometric function, and r = eV ^~ to "> j s a dimensionless scaling parameter so that 
t 3> 1) represents a large (short) time elapsed since the time to. This expression is the first non-trivial result 
of this work. We have obtained an analytical expression for the backscattered current taking into account the effect 
of the nonadiabatic switching of the barrier. In the long time regime this current approaches the stationary value 
I ba (t, -co). 

In order to perform a quantitative analysis of the transient process which is now accessible due to the introduction 
of an abrupt triggering mechanism, we found useful to evaluate the relative change between currents turned on at 
times to and -co, /(r) = j~f|~~~) — !• This is a damped oscillatory function of r with period 2ir and relative 

maxima in t = (2n + 1)tt with n natural. As expected, /(t) goes to zero when r — > co, i.e, when the current is the 
one corresponding to a simple static impurity acting at all times. At this point, in order to have a more intuitive 
picture of the transient process, we define a reference value of r, tr that enables us to identify time scales for which 
the transitory stage is relevant or not. To this end it seems natural to examine the way in which the relative maxima 
of /(t) decrease as r increases. Since the first maximum is located at r = 7r, it is a very good approximation to use 
the following asymptotic (r 3> 1) expression for /(r): 

/(t) w A(t,K) cos[t + 7t], (14) 

where A{t, K) = cos [ 7r ^-] r [ 1 _ 2 jf] i s the envelope of a damped harmonic oscillation. Since the change in the backscattered 
current due to the sudden switching behaves as t~ 2K , we conclude that the relaxation of the system is faster for small 
electron-electron interaction (For the sake of clarity, let us stress that this relaxation characterizes the transition of 
the backscattering current between two off-equilibrium regimes). We thus find an explicit connection between electron 
interactions and the switching time of the externally controlled barrier: the stronger the correlations, the longer the 
persistence of the non adiabatic effect. Formula ([M)) provides a direct way of defining a dimensionless relaxation 
parameter Tr, such that for tr < 1(> 1) the switching process is relevant (irrelevant). We define tr as the value of r 
such that A(t,K) = 1/r (r > 1), which means that for r > tr the value of Ib a (t,to) differs in less than 100/r percent 
from the current obtained when the switching process is not taken into account. To be definite, in the following we 
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FIG. 3: Dimensionless relaxation parameter tr as function of K and for different values of b. 



set r = 10. We then get 

TR= \^K\T\\-2K\) ■ (15) 

The solid line in Figure [3] shows the behavior of tr as function of K at zero temperature. We observe that tr grows 
as lO 1 /^ 2 ^ for strong electron-electron interactions (K — > 0). For weak interactions (K — > 1) the effect of the sudden 
switching is negligible, tr goes to zero as (1 — K) 1 / 2 . 

At this point it is useful to recall that the measurable time in this process is the difference t — to. The time interval 
that it takes the backscattered current to reach a value that differs in exactly 10 percent from the steady current is 
t R — t — ^r-. Using the results shown in Fig. 1 for the zero temperature case, which relate values of K with values 
of tr, we can estimate tR — to- Bias voltages usually applied in experiments with nanodevices range from fiV' s to 
mV's. On the other hand, typical interactions take values in the range .25 < K < .752. With these data wc obtain 
10~ n s < tR — to < 10~ 7 s, which are time intervals that could be measured with recently developed pump-probe 
techniques which reach femtosecond-attosecond time resolutions^ 7 -. 

The relationship between the relaxation process and the strength of Coulombian electron correlations revealed in 
our analysis might provide an alternative way to determine the Luttinger parameter K through time measurements. 
The total current after the switching as a function of time is I(t) — I(oo) + Ct~ 2K cos[J7 1], where il = eV/H, C 
is a constant and I(oo) is the stationary current for a static barrier. A determination of the current as function of 
time with a temporal resolution smaller than 2^ is a direct method to obtain the exponent of the temporal decay, 
and then, the K value of the quantum wire. We emphasize that this proposed method is performed at constant 
source-drain voltage. 

In Figure we illustrate the experimental determination of K according to the ideas mentioned above. The solid 
line represents the measurement of the total current I(t) after the switching of the impurity. All parameters involved 
in this plot (the voltage and the strength of the barrier) are within the range of experimental reach. The temporal 
resolution corresponding to this hypothetical measurement is of the order of the picosecond, which is accessible with 
the already mentioned state of the art pump-probe techniques. The next step is to determine the envelope of the 
curve, given by the dashed line in the figure. This envelope can be fitted with the t~ 2K law in order to find K. 
This procedure, which we depict here for T — can be extended for finite temperatures, as we will show in the next 
section. The power law t~~ 2K is a very good approximation for T <Si . For source-drain voltages of the order of 
mV's, which we used in the case of Figured! this corresponds to T <C lO ^. 
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IV. RESULTS AT FINITE TEMPERATURE 



Taking into account that experiments are performed at very low but finite temperatures, in this section we show 
how our results are affected by thermal effects. Using expression (|11[) for vacuum expectation values of exponentials 
of bosonic fields, we can extend the results of Section III in a straightforward way. The backscattered current with 
the impurity switched on at time — oo (static impurity) becomes: 



M*, -oo) = J^ KT[m (^T,nf^T { K - ^T[K + |] sinh[|-], (16) 

where b = is a dimensionless parameter that characterizes the scale regime in temperature of the system: 6^1 
(<C 1) is the high (low)-temperature regime, for fixed voltage. The perturbative expansion is valid for low temperature 
when ^{^-) K ~ l « 1 and for high temperature when M.(AMl)K-i < 1. 

The backscattered current when the impurity is switched on at a finite time to is: 




f i9(T)r[2A-]sm[rf](B[exp[2H,g,l-2«]- , , 
7 " (M0, = l m- m ±- K] ^ >'^--^ < 17) 

where B is the Incomplete Beta function and * indicates the conjugate of the precedent term. This is the analytical 
generalization of (|13[) at finite temperature. As in the former case, the relative change between currents turned on at 
times to and — oo is a damped oscillatory function of r with the same period and position of relative maxima as before. 
Using the asymptotic approximation, we obtain the following expression for the envelope of the damped oscillation: 

. 4bT[2K]sm[TrK] 
(T ' ' } * (2smh[br])^T[K - ±]T[K + J] sinhlf] ' (1 } 

Thus, the change in the backscattered current due to the sudden switching at finite temperature behaves as 
sinh[&T] _2K = sinh[7rfcf,r(i — to)/ft\~ 2K ■ For high temperatures, (b ^> 1), the decay is faster and becomes expo- 

nential: A(r, K, b) a g^ygl^l . 

Experimentally, the high temperature regime can be accessed for low voltages. Temperatures involved in the 
crossover regime (b ~ 1) are T ~ 10° K and T ~ 0.01°A", for applied voltages of order lmU and lfiV, respectively. 
In general, if we know the temperature of the system, a method completely analogous to the one described at the 
end of the previous section can be established to obtain the parameter K . The total current after the switching is 
I(t) = /(oo) + C swh[irkbT(t — to)/h]~ 2K cos[eV t/H]. Thus, the value of K can be obtained by measuring the current 
as a function of time at fixed voltage. 

As in the case of zero temperature, we can define a dimensionless relaxation parameter tr in the same fashion; 
i.e. by determining the value of r such that A(t, K, b) =0.1. We obtain for tr the following general expression as 
function of temperature and electron-electron interaction strength: 



A0bT[2K] sin[7r#] \ 
T[K-4 E }T[K+4 E ]su&[%]) \ ■ (19) 

Figure [3] shows the dimensionless relaxation parameter tr of the system as function of K, for different values of 
temperature. We observe that the distortion in the backscattered current caused by the transient process is more 
important for high electron-electron interactions (K <C 1) and in the regime of zero or low temperature (b <C 1), 
where tr^> 1. In this case the relaxation is very slow and the sudden switching changes significantly the value of the 
backscattered current for a long time. If b and K grow, this effect tends to disappear. For high temperatures, the 
switching effect becomes irrelevant. 



tr = — arcsinh 
2b 



V. CONCLUSIONS 



We have theoretically analyzed the consequences of a sudden switch of a tunnel barrier on transport properties, 
in a long and clean LL. In particular we studied the behavior of the backscattered current h s in a system subjected 
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to a bias voltage V. The backscattering time-dependent impurity is assumed to be point-like and weak. Under these 
conditions, using bosonization and our result (jlip for the expectation value of exponentials of bosonic operators, we 
obtained an analytical expression for Ib s as a function of r = eV ^ to ^> (f Q [ s the instant at which the barrier is switched 
on), b = if and the Luttinger parameter K (See formula (fTTj) ). At long times t ^> to, hs(t,to) reaches a steady- 
state value that coincides with the value of the current corresponding to a static impurity (a barrier switched on at 
an infinitely remote instant in the past) l(, s (i, — oo). By carefully examining the way Ib s (t,to) approaches Ib s (t, — oo) 
as time goes by, we intuitively characterized the transient process defining a dimensionless relaxation parameter tr. 
Employing an asymptotic expression for the exact solution (JTTJ) , we obtained a simple expression for tr as function of 
b and K (See formula (fl"9f and Fig. 2). From this result one concludes that, for fixed bias, tr grows with decreasing 
K, meaning that transient effects on the current are expected to last longer when the electron correlation is higher. 
In fact, tr — > for a free system (K — > 1). One can also see that the transitory effect is more pronounced for low 
temperatures (b <C 1). 

Finally, we would like to stress that our study could have some interesting experimental applications. Indeed, from 
our main result (|17[) , and taking into account that the total current after the switching is I{t) = I — I{, s (t,tQ), we 
obtained I(t) = 7(oo) + C sinh[7rfcf,T(t — to)/h]~ 2K cos[eVt/h}. This means that, at fixed temperature and bias voltage, 
the value of K of a ID structure could be determined by performing measurements of the total current as function 
of time, after turning on a tunnel barrier through a localized gate voltage, following, for instance, the techniques of 
reference^. 
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